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We present measurements and theoretical interpretation of the magnetic field depen- 
dent excitation spectra of a two-electron quantum dot. The quantum dot is based on 
an ALjGai— x As parabolic quantum well with effective g*-factor close to zero. Results 
of tunneling spectroscopy of the four lowest states are compared to exact diagonaliza- 
tion calculations and a generalized Heitler-London approximation and good agreement 
is found. Electronic correlations, associated with the formation of an H2-type Wigner 
molecule, turn out to be very important in this system. 
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1. Introduction 

The attempt to understand the effects of Coulomb interaction in solids — with 
exchange- and correlation effects being the most interesting manifestations — has 
always been a driving force behind research in solid state systems. Nowadays, fabri- 
cation techniques for quantum dots offer the unique possibility to tailor interacting 
quantum systems to an unprecedented degree. Understanding quantum dot helium, 
a man-made two-electron system, is a hallmark for our ability to design and control 
more complex interacting quantum systems such as those required for the imple- 
mentation of quantum information processing schemes. Here, we present a study of 
quantum dot Helium fabricated on a parabolic quantum well 1 which goes beyond 
earlier work 2 by detecting higher lying excited states at high magnetic fields and 
by presenting evidence for the importance of correlation effects. 



I 
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2. Parabolic Quantum Wells 

Parabolic quantum wells (PQWs) based on GaAs/Al x Gai_ x As heterostructures 
have been introduced in 1988 for creating high mobility quasi-three-dimensional 
electron gases. 3 Within the effective mass approximation their conduction and va- 
lence band edges vary parabolically along the growth direction, as schematically 
depicted in Fig. 1(a), caused by the varying Al content x of the Al x Gai_ x As mate- 
rial. For example, the conduction band edge is described by 



where a is the curvature parameter of the parabolic confinement. 

A PQW can be filled with electrons using standard remote doping techniques. 
The electrons occupy subbands which form in the well self-consistently under the in- 
fluence of the parabolic confinement and screening. With increasing two-dimensional 
electron density, this interplay leads — in the limit of large density — to a constant 
three-dimensional electron density — eeaa/e along with an increasing width of 
the density distribution in growth direction. Here, e is the relative dielectric constant 
of the material. 

A very peculiar property of PQWs is the rigidity of the subband wave functions 
under the influence of an external homogeneous electric field E z in growth direction, 
given that the two-dimensional electron density is constant. Such a field can, for 
example, be created with metallic gates on top and at the back side of the substrate. 
It is equivalent to an electrostatic potential changing linearly with z which adds to 
the built-in z 2 -confinement in Eq. (1) and leads to a net parabolic potential with 
the same curvature parameter a, but shifted by Az = eE z /a in growth direction. If 
this shifted potential is populated with the same two-dimensional electron density, 
it will develop the same self-consistent subbands with the same envelope functions, 
but shifted in space by Az. This property of PQWs has been exploited for mapping 
the subband wave function in real space. 4 

In the Al x Gai_ x As material, the effective mass m* and the effective g-factor g* 
depend on the Al fraction x. For example, the effective g-factor varies according to 



The wave functions of subbands in PQWs probe regions of varying x and therefore 
feel an average g*. This average varies, when the wave function is shifted, providing 
a gate voltage tunable g* , i.e., a gate voltage tunable Zeeman splitting in a magnetic 
field. 5 Furthermore, PQWs can be designed to have certain g* values when the wave 
function is positioned symmetrically around the minimum of the parabola at zero 
gate voltages by deliberately choosing a certain offset value for x at the well center. 
For example, PQWs with g* ss can be grown. 

Recently, PQWs have served as the basic material for the fabrication of meso- 
scopic devices such as quantum point contacts 6 and quantum dots. 1,7 ' 8 In particu- 
lar, the spin states in a quantum dot fabricated on a PQW were investigated 8 and 




(1) 



g * = -0.44 + 3.833a; for < x < 0.3. 



(2) 
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few-electron dots with 0-3 electrons were realized. 1 Such few-electron dots are can- 
didates for the implementation of spin-qubits for quantum information processing. 9 
Dots on PQWs promise the additional feature of a tunable g* which could be used 
for the realization of gate-controlled single-qubit operations. 10 

3. Sample and Experimental Setup 

The sample used for this study 1 is a quantum dot (QD) with integrated charge 
readout fabricated on a PQW by electron-beam lithography. The PQW material 
has been designed and measured to have g* ~ 0. A highly doped layer, 1.3 //m 
below the well serves as a back gate. A scanning electron micrograph of the sample 
is shown in Fig. 1(b). The QD is formed by applying negative voltages to the gates 
PCI, PC2, PG and MG. The quantum point contact formed between gates PC2 
and CD has been used as a detector for the charge in the quantum dot 11 and 
allows to count the number of electrons in the dot from upwards. Experiments 
were performed in a dilution refrigerator with a base temperature of 100 mK using 
standard low-frequency DC and AC conductance measurement techniques. 

4. Excited State Spectroscopy of Quantum Dot Helium 

The differential conductance of the QD has been measured in the few-electron 
regime as a function of plunger gate voltage V pg and source-drain voltage Vbias- 
Figure 2 shows the resulting Coulomb-blockade diamonds with the electron numbers 
N indicated. The QD has a charging energy Afi 2 = 6.9 meV for the second electron. 
The lowest excited state ES of the 1-electron system has an energy Ai = 5 meV 
above the ground state. Therefore, confinement and interaction effects are of com- 
parable order of magnitude in this system. Singlet (S) and Triplet (T) resonances 
are observed in the differential conductance of the 2-electron dot. The singlet-triplet 



MG 




Fig. 1. (a) Schematic plot of the conduction band edge in a parabolic quantum well as a function 
of coordinate z measured in the direction of crystal growth. The lowest bound subband state 
with its envelope wave function is indicated, (b) Scanning electron micrograph of the few-electron 
quantum dot with integrated charge readout fabricated on a PQW. 
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v bias ( mV ) magnetic field (T) 



Fig. 2. (a) Differential conductance of the few-electron quantum dot as a function of the plunger 
gate voltage Vp g and the source-drain voltage Vbi as - Electron numbers are indicated. ES labels the 
resonance of the first excited state of the 1-electron dot, S and T label the singlet- and triplet state 
transitions of the 2-electron system, (b) The top panel shows the measured excitation spectrum 
of QD He in a magnetic field. The bottom panel is the corresponding energy spectrum of QD He 
calculated with EXD. The gray-shaded region indicates the bias- window of Vbi as = 2.5 meV of the 
experiment. 

separation J is about 2 meV. 

For the investigation of the excitation spectrum of quantum dot Helium, i.e., 
the two-electron system, in a magnetic field, the QD was tuned into a slightly more 
open regime allowing the plunger gate to cover the range from N = 1 to 3. The 
excitation spectrum of the system was then recorded at fixed source-drain voltage 
Vbias = 2.5 mV as a function of plunger gate voltage V pg and magnetic field B. The 
result 1 is shown in Fig. 2(b) (upper panel). A number of transitions can be followed 
as a function of magnetic field: most prominent is the transition from the 1-electron 
ground state to the 2-electron singlet state, labeled So- The first excitation at low 
B is the triplet state labeled T + which becomes the ground state at a magnetic field 
B ss 4 T. The transition labeled S2 appears at B w 3 T in the bias window and 
shows an avoided crossing with Sq indicating its spin-singlet character. The state 
T +! cm appearing at even higher B can only be identified by comparing to the results 
of an exact diagonalization calculation (EXD) to be described below. It turns out 
to be combined excitation of the center of mass (CM) motion, the relative motion 
(+) and the spin degree-of-freedom (T). 

5. Theoretical Methods and Calculations for Quantum Dot Helium 

In order to interpret the measured excitation spectra in detail, we present an exact 
diagonalization (EXD) and an approximate (generalized Heitler-London, GHL) mi- 
croscopic treatment for two electrons in a single elliptic QD specified by parameters 
that correspond to our experimental device. 1 
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The Hamiltonian for the two 2D interacting electrons is given by 

U = H{v 1 ) + H{v 2 )+ 1 e 2 /{Kr 12 ), (3) 

where the last term is the Coulomb repulsion, n (12.5 for GaAs) is the dielectric 
constant, and r\ 2 = \v\ — r 2 |. The prefactor 7 accounts for the reduction of the 
Coulomb strength due to the finite thickness of the electron layer in the z direction 
and for any additional screening effects due to the gate electrons. H (r) is the single- 
particle Hamiltonian for an electron in an external perpendicular magnetic field B 
and an appropriate confinement potential. For an elliptic QD, the single-particle 
Hamiltonian is written as 

H(r)=T+±m*(Lu 2 x x 2 +u J y), (4) 

where T = (p - eA/c) 2 /2m*, with A = 0.5(— By, Bx, 0) being the vector potential 
in the symmetric gauge. The effective mass is m* = 0.07too, and p is the linear 
momentum of the electron. The second term is the external confining potential. In 
the Hamiltonian (4), we neglect the Zeeman contribution due to the negligible value 
(g* « 0) of the effective Lande factor in our sample. 

5.1. Generalized Heitler London Approach 

The GHL method for solving the Hamiltoninian (3) consists of two steps. In the 
first step, we solve self-consistently the ensuing unrestricted Hartree-Fock (UHF) 
equations allowing for lifting of the double-occupancy requirement (imposing this 
requirement gives the restricted HF method, RHF). For the S z = solution, this 
step produces two single-electron orbitals u L ^ R (r) that are localized left (L) and 
right (R) of the center of the QD [unlike the RHF method that gives a single 
doubly-occupied elliptic (and symmetric about the origin) orbital]. At this step, the 
many-body wave function is a single Slater determinant v I'uhf(1 T,2 J.) = \ur,(l f 
)ur(2 J.)) made out of the two occupied UHF spin-orbitals Ul(1 |) = wi,(ri)a(l) 
and u R (2 J.) = u R (r 2 )/3(2), where a{(3) denotes the up (down) [f (J,)] spin. This UHF 
determinant is an eigenfunction of the projection S z of the total spin S = §1 + §2, 
but not of S 2 (or the parity space-reflection operator). 

In the second step, we restore the broken parity and total-spin symmetries by 
applying to the UHF determinant the projection operator 12 ' 15 Vgp in — 1 T ^12, 
where the operator zui 2 interchanges the spins of the two electrons; the upper (mi- 
nus) sign corresponds to the singlet. The final result is a generalized Heitler-London 
two-electron wave function ^Q HL (ri,r 2 ) for the ground-state singlet (index s) and 
first-excited triplet (index t), which uses the UHF localized orbitals, 

^Hi>i,r2) oc (u L (r 1 )u R (r 2 )±u L (r 2 )u R {r 1 )) X s ' t , (5) 

where x s '* = (a(l)/3(2) =Fa(2)/3(l)) is the spin function for the 2e singlet and triplet 
states. The general formalism of the 2D UHF equations and of the subsequent 
restoration of broken spin symmetries can be found in Rcfs. 12 > 13 > 14 > 15 . 



4, 2008 13:51 WSPC/INSTRUCTION FILE Ihn 



6 T. Ihn, C. Ellenberger, C. Yannouleas, U. Landman, K. Ensslin, D. Driscoll, andA.C. Gossard 

The use of optimized UHF orbitals in the GHL is suitable for treating sin- 
gle elongated QDs. The GHL is equally applicable to double QDs with arbitrary 
interdot-tunneling coupling. 12,15 In contrast, the Heitler-London (HL) treatment 16 
(known also as Valence bond) , where non-optimized "atomic" orbitals of two iso- 
lated QDs are used, is appropriate only for the case of a double dot with small 
interdot-tunneling coupling. 17 

The orbitals Ul,r{t) are expanded in a real Cartesian harmonic-oscillator basis, 

i.e., 

ML , fl (r)=]T^V,(r), (6) 

where the index j = (m, n) and <fj(r) = X m (x)Y n (y), with X m (Y n ) being the 
cigenfunctions of the one-dimensional oscillator in the x(y) direction with frequency 
u> x (u)y). The parity operator V yields VX m (x) = {— l) m X m {x), and similarly for 
Y n (y). The expansion coefficients C. ' are real for B — and complex for finite B. 
In the calculations we use K = 54 and/or K = 79, yielding convergent results. 

5.2. Exact Diagonalization 

In the EXD method, the many-body wave function is written as a linear superpo- 
sition over the basis of non- interacting two-electron determinants, i.e., 

IK 

^ D (r 1 ,r 2 ) = ^^/|^(l;iM2;j)), (7) 

i<j 

where tp(l; i) = <^(1 |) if 1 < i < K and ip(l;i) = ipi- K {l I) if K + 1 < i < 2K 
[and similarly for ip(2,j)]. The total energies -EgXD ano - tne coefficients f2*j are 
obtained through a "brute force" diagonalization of the matrix eigenvalue equation 
corresponding to the Hamiltonian in Eq. (3). The EXD wave function does not 
immediately reveal any particular form, although, our calculations below show that 
it can be approximated by a GHL wave function in the case of the elliptic dot under 
consideration. 

5.3. Results and Comparison with Measurements 

To model the experimental quantum dot device, we take, following Rcf. 1, hw x = 
4.23 meV, hw y — 5.84 meV, and 7 = 0.862. The corresponding anisotropy is 
ujy/uj x — 1.38, indicating that the quantum dot considered here is closer to be- 
ing circular than in other experimental systems. 18,19 

As shown in Rcf. 1 and Fig. 2(b), the experimental findings can be quantitatively 
interpreted by comparing to the results of the EXD calculations for two electrons 
in an anisotropic harmonic confinement potential with the parameters listed above. 
All the states observed in the measured spectra (as a function of the magnetic field) 
can be unambiguously identified 1 with calculated ground-state and excited states of 
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Fig. 3. (a) Comparison of J(B) calculated with different methods and the experimental results 
(open squares). Solid line: EXD. Dotted line: GHL. Dashed line: RHF. (b) Single-particle UHF 
orbitals (modulus square) that are used in the construction of the GHL wave function in Eq. (5). 
Lengths in nm and orbital densities in 10 -3 mi" 2 . Arrows indicate up and down spins. For the 
parameters used in the calculation to model the anisotropic QD, see text. 

the two-electron Hamiltonian [see Fig. 2(b)]. In particular, the calculated magnetic- 
field-dependent energy splitting, Jexd(-S) = E^ XD (B) — E^ XD (B), between the two 
lowest singlet (.So) and triplet (T + ) states is found to be in remarkable agreement 
with the experiment [see Fig. 3(a)]. 

A deeper understanding of the structure of the many-body wave function can 
be acquired by comparing the measured J(B) with that calculated within the GHL 
and RHF approximations. To facilitate the comparisons, the calculated Jghl(B) 
and Jrhf(-B) curves are plotted also in Fig. 3(a), along with the EXD result and the 
measurements. Both the RHF and GHL schemes are appealing intuitively, because 
they minimize the total energy using single-particle orbitals. It is evident, however, 
from Fig. 3(a) that the RHF method, which assumes that both electrons occupy a 
common single-particle orbital, is not able to reproduce the experimental findings. 
On the contrary, the GHL approach, which allows the two electrons to occupy two 
spatially separated orbitals, appears to be a good approximation. Plotting the two 
GHL orbitals [see Fig. 3(b)] for the singlet state clearly demonstrates that the two 
electrons do not occupy the same orbital, but rather fill states that are spatially 
separated significantly. 

The UHF orbitals from which the GHL singlet state is constructed [see Eq. (5)] 
are displayed on Fig. 3(b) for both the B = and B = 3.8 T cases. The spatial 
shrinking of these orbitals at the higher £?-value illustrates the "dissociation" of the 
electron dimcr with increasing magnetic field. The asymptotic convergence (beyond 
the ST point) of the energies of the singlet and triplet states, [i.e., J(B) — > as 
B — > oo] is a reflection of the dissociation of the 2e molecule, since the ground-state 
energy of two fully spatially separated electrons (zero overlap) does not depend 
on the total spin. We stress again that the RHF, which corresponds to the more 
familiar physical picture of a QD-Hclium atom, fails to describe this dissociation, 




(b) 



B=0 
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B=0 CPDs B=3.8 T 




Fig. 4. CPDs extracted from the EXD wave function for the singlet state for B = and B = 3.8 
T. The CPD expresses the conditional probability for finding the second electron at position r 
given that the first electron is located at ro (denoted by a heavy solid dot). For the parameters 
used in the calculation to model the anisotropic QD, see text. Lengths in nm and CPDs in arbitrary 
units. 

because Jrhf (B) diverges as the value of the magnetic held increases. 

In contrast to the RHF, the GHL wave function is able to capture the impor- 
tance of correlation effects. Further insight into the importance of correlations in 
our QD device can be gained through inspection 1 of the conditional probability 
distributions 20 (CPDs) associated with the EXD solutions; see an illustration in 
Fig. 4. Indeed, already at zero magnetic held, the calculated CPDs provide further 
support of the physical picture of two localized electrons forming a state resembling 
an H^-type 12 ' 1 Wigner molecule. 21 ' 22 

5.4. Degree of Entanglement 

Further connections between the strong correlations found in our microscopic treat- 
ment and the theory of quantum computing 17 can be made through specihcation 
of the degree of entanglement between the two localized electrons in the molecular 
dimer. For two electrons, we can quantify the degree of entanglement by calculating 
a well-known measure of entanglement such as the von Neumann entropy 23 ' 2 for 
indistinguishable particles. To this end, one needs to bring the EXD wave function 
into a diagonal form (the socalled "canonical form" 23 ' 25 ), i.e., 

M 

vf^ D (n,r 2 ) = Y, 4' l \m 2fc - 1)*(2; 2k)), (8) 

k=l 

with the <j>(i)'s being appropriate spin orbitals resulting from a unitary transforma- 
tion of the basis spin orbitals ^(j')'s [see Eq. (7)]; only terms with Zk ^ contribute. 
The upper bound M can be smaller (but not larger) than K (the dimension of the 
single-particle basis); M is referred to as the Slater rank. One obtains the coeffi- 
cients of the canonical expansion from the fact that the \zk\ 2 are eigenvalues of the 
hermitian matrix Q^Q [fi, see Eq. (7), is antisymmetric]. The von Neumann entropy 
is given by S — — X)fc=i l z fc| 2 log 2 ( | | 2 ) with the normalization 52k=i \ z k\ 2 = 1< 

The EXD singlet has obviously a Slater rank M > 2. The von Neumann en- 
tropy for the EXD singlet (5|- XD ) is displayed in Fig. 5. It is remarkable that 
iSexd increases with increasing B, but remains close to unity for large B, although 
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Fig. 5. Von Neumann entropy for the lowest singlet EXD state of the elliptic dot as a function 
of the magnetic field B. On the top, we show histograms for the \zk\ 2 coefficients [see Eq. (8)] 
of the singlet state at B = 3.8 T (left) and B = 8.0 T (right) illustrating the dominance of two 
determinantal configurations (in agreement with the GHL picture). Note the small third coefficient 
| ^3 1 2 = 0.081 for B = 8.0 T. For the parameters used to model our device, see text. 

the maximum allowed mathematical value is \og 2 (K) [for example for K = 79, 
log 2 (79) = 6.3]. The saturation of the entropy for large B to a value close to unity 
reflects the dominant (and roughly equal at large B) weight of two configurations in 
the canonical expansion [see Eq. (8)] of the EXD wave function, which are related 24 
to the two terms in the canonical expansion of the GHL singlet. This is illustrated 
by the histograms of the \zf.\ 2 coefficients for B = 3.8 T and B = 8.0 T in Fig. 5 
(top). Notice that the ratio \z 2 \ 2 /\zi\ 2 reflects the extent of the overlap between the 
two GHL orbitals 24 , with the ratio increasing for smaller overlaps (corresponding 
to a more complete dissociation of the Wigner molecule). 

The above discussion illustrates that microscopic calculations that are shown to 
reproduce experimental spectra 1 can be used to extract valuable information that 
allows assessment of the suitability of a given device for quantum computations. 



6. Conclusions 

Measurements and theoretical interpretation were presented of the magnetic field 
dependent excitation spectra of a two-electron quantum dot. The quantum dot is 
based on an ALGai-^As parabolic quantum well with effective g*-factor close to 
zero. Results of tunneling spectroscopy of the four lowest states were compared 
to exact diagonalization calculations and a generalized Hcitlcr -London approxima- 
tion and good agreement was found. Electronic correlations, associated with the 
formation of an H 2 -type Wigner molecule, turn out to be very important in this 
system. 
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